Nontrivial interplay between superconductivity and spin-orbit coupling in non-centrosymmetric 

ferromagnets 
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Motivated by the recent discoveries of ferromagnetic and non-centrosymmetric superconductors, we present 
a mean-field theory for a superconductor that both lacks inversion symmetry and displays ferromagnetism, a 
scenario which is believed to be realized in Ulr. We study the interplay between the order parameters to clarify 
how superconductivity is affected by the presence of ferromagnetism and spin-orbit coupling. One of our key 
findings is that the spin-orbit coupling seems to enhance both ferromagnetism and superconductivity in all spin 
channels. We discuss our results in the context of the heavy fermion superconductor Ulr and analyze possible 
symmetries of the order parameter by the group theory method. 

PACS numbers: 74.20.Rp, 74.50.+r, 74.20.-z 
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In the past decade, a number of superconductors have been 
discovered that are called 'unconventional' as they fall outside 
the Bardeen-Cooper-Schrieffer (BCS) paradigm of electron- 
phonon mediated pairing with an isotropic gap. Of those, 
UPt3 ' and Sr2Ru04^ were the first compounds to have been 
confirmed as triplet p-wave superconductors. More recently, a 
ferromagnetic (FM) superconductor was discovered in UGe2 
under pressure^, where the presence of an internal FM mo- 
ment strongly suggests that only the equal-spin triplet pair- 
ing survives. In this latter example both the time-reversal and 
the gauge symmetry due to SC order are spontaneously bro- 
ken, which made UGe2, as well as its cousins URhGe'* and 
UCoGe^ an exciting avenue for theoretical and experimental 
research. 

For spin-triplet pairing, Anderson noticed^ that inversion 
symmetry is required to obtain a pair of degenerate states 
cj^|0) and c^j^jO) capable of forming a Cooper pair It was 
therefore surprising that superconductivity was discovered in 
the heavy fermion compound CePtaSi which lacks inversion 
symmetry It soon became clear however that in the case of a 
non-centrosymmetric crystal, the spin-orbit coupling (SOC) 
mixes different spin states, so that the division into triplet 
and singlet symmetry of the SC order parameter becomes 
meaningless. A bulk of theoretical work exists that has pro- 
vided a symmetry -based phenomenology to explain this in de- 
tail a*^'^'^"'"'' ^. The symmetry of the superconducting (SC) gap 
in this and other unconventional superconductors is presently 
a matter of intense investigatio n'^-'"^''^-'^ . 

An intriguing question is what happens if time-reversal 
symmetry is broken in a crystal that lacks a centre of inver- 
sion. Can such a material become a superconductor? This 
question was answered affirmatively when superconductiv- 
ity was discovered in the non-centrosymmetic ferromagnetic 
compound Ulr under pressure'^. The symmetry of the SC or- 
der parameter and its connection to FM nevertheless remains 
unclear, which motivates the present study. Spontaneous sym- 
metry breaking in condensed matter systems is conceptually 
of immense importance, as it may provide clues for what 
could be expected in systems belonging to vastly different ar- 
eas of physics. The study of a condensed-matter system such 



as Ulr with multiple broken symmetries is likely to have im- 
pact on a number of disciplines of physics, including such 
disparate phenomena as mass differences between elementary 
particles and extremely dilute ultra-cold atomic gases. 

In this work, we study a model system of a non- 
centrosymmetric superconductor with substantial spin-orbit 
coupling, which at the same time exhibits itinerant ferromag- 
netism. The origin of the SOC may be either that the crystal 
structure lacks a center of inversion, such as in Ulr, or due to a 
thin-film geometry where the breakdown of inversion symme- 
try near the surface induces transverse electrical fields, leading 
to the well-known Rashba SOC Our model should therefore 
be relevant both to the non-centrosymmetric and centrosym- 
metric heavy fermion compounds, since the SOC is consider- 
able in any case due to the high atomic number Specifically, 
materials that exhibit coexistence of SC and FM order and 
where SOC is large include UGejl, URhGel, UCoGe?, and 
Ulr—. For this model, we construct a mean-field theory, solve 
the saddle point equations for the order parameters and study 
the effect of spin-orbit coupling on the superconducting order 
parameters. Finally, we discuss application of this model to 
the case of Ulr. 

To label the SOCh-FM spUt bands, it is possible to introduce 
a pseudospin basis in which the normal-state Hamiltonian is 
diagonalized. In the original spin basis, the SC matrix order 
parameter is characterized, in analogy to the p-wave stato^^, 
by a vector dk and scalar As so that AQ^(k) = iA^CTj, + 
[i(dk • cr)ay\ap. Note that, unlike the usual p-wave SC, a 
singlet component Ag of the gap will also be present since 
antisymmetric SOC in general mixes the parity of the order 
parameter. Below, . ? . is used for 2x2 matrices. We now 
proceed to write down the effective Hamiltonian H = + 
Hsc for our system. In the normal state, the Hamiltonian in 
momentum-space reads^^ 

^ Hq + [CkalSki - hd-^+a ■ gk)Q/3Ck/3], (1) 
kQ/3 

where Hq = /iVM^/2. Above, the dispersion relation £k is 
measured from chemical potential p, and the magnetization 
M — |M| is taken along the easy-axis, while h = lAI is 
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the exchange splitting of the bands and gk is the SOC vector. 
When superconductivity coexists with FM, the SC pairing is 
generally believed to be non-unitary22,, characterized by dk x 
d^. 7^ 0. In such a scenario, the SC order parameter couples to 
the spontaneous magnetization M through a term 7M • dk x 
dj^ in the free energy, where the sign of 7 is determined by the 
gradient of the DOS at Fermi level^i and (Sk) = idk x dj^ 
is the spin associated with the Cooper pair Thus, for 7 < 
it is expected that a SC pairing state obeying idk x dj. || M 
is energetically favored, implying that dk must be complex- 
valued. Our model captures broken time-reversal symmetry in 
addition to antisymmetric SOC. As shown by Anderson^, the 
presence of the latter is detrimental to spin-triplet SC pairing 
state, unless dk || gk- In our case, it is obvious that a non- 
unitary SC pairing state cannot satisfy this condition since dk 
is complex, whereas gk must be real for the Hamiltonian to 
be hermitian. 

The SOC vector reads gk = — g-k, and we introduce 
ffk = Sk.x — igk.a for later use. We consider the SOC in 
the Rashba form, namely gk = X{ky,—kx,0)- This cor- 
responds to a situation where an asymmetric potential gradient 
is present along the z-axis, and is also the scenario realized in 
non-centrosymmetric CePtsSi" . We have introduced fermion 
operators {ckcr} in a basis ip]^ = [0^1,0^1]'^. 

Diagonalizing the normal-state Hamiltonian yields the 
quasiparticle excitations i^kcr = — a-y/h^ + A^fc^, which 
due to the SOC are characterized by the pseudospin a = 
±1. For later use, we define A/k = [1 + \'^k^/{h + 
Y^7^2~j7^^2^-j2j-i/2 rpj^g superconducting pairing is now as- 
sumed to occur between the excitations described by (^k- Due 
to the presence of antisymmetric spin-orbit coupling, this au- 
tomatically leads to a mixed-parity SC state in the original 
spin basis. To see this, we introduce 
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and perform a standard mean-field decoupling, which after 
an additional diagonalization yields the total Hamiltonian in 
the superconducting state: H — Hq + Eko-(£'ko- — E^a — 
Aka^L + 2ryL'7k.)/2, where i?k. = (^L + lAk.P)^/' 
and {?7kCT, ijka} are new fermion operators. The merit of this 
procedure is that we can now obtain simple self-consistency 
equations for the gaps Akcr, which may then be transformed 
back to the gaps in the original spin-basis (^k by means of 
the unitary transformation Pk- We assume a chiral p-wave 
symmetry for the gaps with a corresponding pairing potential 
Vkk'cr = — ffscC"^^'^^'^ ^ where tan = k^/ky. The moti- 
vation for this is that this choice ensures that the condition 
dk II gk is satisfied exactly for h ^ Q, and corresponds to 
a fully gapped Fermi surface which favors the condensation 
energy. The gaps obtain the form Aku = — crAcr^oe"^'^ and 
we find a self-consistency equation of the standard BCS form 
with a cutoff uj on the pairing-fluctuation spectrum which we 
do not specify further Moreover, N'^ {e) is the pseudospin- 
resolved density of states (DOS) for the E'kcr (o'=±) bands 
of the quasiparticle excitations^'*. Introducing the total DOS 
at the Fermi level for a normal metal Nq = mV\/2mii/TT^ 




FIG. 1: (Color online) Self-consistent solution of the order param- 
eters (a-d) as a function of the FM exchange parameter /, the ratio 
between the singlet and triplet gaps (e) _Ra ~ A||/(A| + A|) and 
the maximal critical temperature (f) Tc, max as a function of the FM 
exchange parameter /. 



and defining c = gNo/2, the analytical solution for the gaps 
reads A,,o = 2wexp{-l/[ci?„(0)]}, i?,(e) = 2N''{e)/No. 
With the analytical solution for Acr.o in hand, we may exploit 
the unitary transformation Pk to express the superconducting 
gaps in the original spin basis as follows: 



^kT 



^ki 



^kTi 



:'^[AT,o(<)^ + Ai,o(AAk^)'A2fc2(o)A: 



-^A,,o(AA^)'A|fc,(0)|Ak., fT = ±l, (3) 



where we have defined JV-^ = A/'k=k''(o) and Ako- — [h + 
^h? + A^fc^(O)]^"'". Note that in the original spin basis, 
the superconducting order parameter is in general a mixture 
of triplet (Akcr) and singlet (Ak|x) components. The self- 
consistency equation for the magnetization is: 
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aAeR''{e)he 



[/i2 + A2fc2(e)](e2+A^^„) 



= 0, (4) 



where the integration is over the bandwidth and / = /Aq. 
Eqs. ([3] nil are the main analytical results of this work. 

Let us briefly investigate some important limiting cases of 
Eq. (|3]l. In the absence of spin-orbit coupling (A 0), one 
finds 1 and Akcr = Akcr while Ak|x — 0, such that 

we reproduce the results of Refsi22i21. In the absence of an 
exchange energy (h 0), one finds that 1/V2 and 

AkT = -e"^(A|^o + Ax,o)/2, Akx - e-'^(A|^o + Ax,o)/2, 
and Ak|x = (Aj^.o — A| o)/2. As demanded by consis- 
tency, the triplet gaps are equal in magnitude since there is 
no exchange field and the singlet component is nonzero since 
A| 7^ Aj^ in general. Finally, Eq. (|4|i reproduces the well- 
known Stoner criterion / > 1 for the onset of FM in the ab- 
sence of SOC and SC (A ^ 0, 5 ^ 0). 

We now focus on the general case in which h and 
A 7^ 0. First of all, we must specify the range of the pa- 
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rameters in the problem that corresponds to a physically re- 
alistic scenario. We allow h to range, in principle, from 
to fi, the latter denoting a fully polarized ferromagnet. As 
a convenient measure of the strength of SOC, we introduce 
the dimensionless quantity a^oc = ■\/2A^m//i which has a 
direct physical interpretation: namely, it is the ratio of the 
SOC (at Ep) to the Fermi energy fi. The parameter a^oc is 
allowed to vary from to 6, where 6 denotes a fraction of 
the Fermi energy. We take (5 = 0.5 as a sensible upper limit. 
Note that generically, the SOC strength at the Fermi level is 
different for the two quasiparticle bands, and moreover de- 
pends on h. For a given value of h, one may derive that 
A < Sfi/[2fim+ y^2m'^{h'^ + 6^ii^)]^^^ ensures that the spin- 
orbit energy is less than S x jj. for both quasiparticle bands. 

In Fig. [iH-d, we present the self-consistent solutions for 
the order parameters in Eqs. (O and (|4| as a function of the 
FM exchange parameter / for several values of asoc- We have 
defined Ao- = lAkal and A|| = |Akti|, and fixed lo/ ji = 
0.01 and m/ jj. = 5 x 10* with c = 0.2, which are standard 
choices. For Ofsoc = 0, the onset of FM occurs at / = 1.0 
which lifts the degeneracy of A| and A|, while A| j is always 
zero. Upon increasing a^oc, it is interesting to note that the 
PM-FM transition occurs at lower values of /, indicating that 
spin-orbit coupling favors ferromagnetic ordering. For a^oc 7^ 
0, it is seen that A|| is also non-zero, although it becomes 
suppressed at the onset of ferromagnetism. A common feature 
for all gaps is that they increase with Ofsoc in the absence of 
ferromagnetism and deep inside the ferromagnetic phase / > 
1.02. In the intermediate regime, there are crossovers between 
the gaps for different values of a^oc due to the different onsets 
of ferromagnetic order. By comparing the behaviour between 
the gaps for increasing / with ajoc 7^ 0, one infers that A| and 
A|j eventually saturate at a constant non-zero value, while 
A I continues to increase steadily. This is quite different from 
the case when a^oc = 0, where the minority spin-gap goes 
to zero rapidly with increasing /. This seems to suggest that 
the presence of spin-orbit coupling in the system ensures the 
survival of the minority-spin gap A| and the singlet gap A|| 
even though the FM exchange energy becomes strong. 

In Fig. [T^ and[T]F, we plot the ratio of the singlet and triplet 
gaps, defined as i?A — AtJ./(^T + ^i)' maximal 
critical temperature Tcmax for the onset of superconductiv- 
ity. It is seen from the left panel that i?A increases with a^oc 
in the PM regime, suggesting that the singlet component be- 
comes more prominent in the system as compared to the triplet 
gaps. However, at the onset of FM order, i?A decreases since 
the singlet component becomes suppressed by the Zeeman- 
splitting. In the right panel, one observes that Tcmax in- 
creases both with ctsoc and /. Our findings suggest that the 
presence of antisymmetric SOC, originating from e.g. non- 
centrosymmetricity of the crystal structure, enhances both the 
tendency towards ferromagnetism and the magnitude of the 
SC gaps in all spin channels. In the absence of spin-orbit cou- 
pling, it was shown in Ref.^^ that the simultaneous coexistence 
of FM and non-unitary triplet superconductivity is the ther- 
modynamically favored state as compared to the pure normal, 
FM, or SC state. Since the presence of spin-orbit coupling is 
seen to enhance both the FM and SC order parameters, it is 



reasonable to expect that the coexistent state is still thermody- 
namically the most favorable one even when asoc 7^ 0. 

Out of the known non-centrosymmetric superconductors, 
Ulr is the only compound that is also a ferromagnet. This ma- 
terial, which is ferromagnetic at ambient pressure, develops 
superconductivity in a narrow pressure region around P ^ 
2.6 GPa right next to the FM-PM quantum phase transition, 
with a maximum SC transition temperature Tsc ^ 0.14 K'^. 
At this pressure, the saturated magnetic moment was mea- 
sured to be 0.07 /is per U atom, and such a small value clearly 
indicates the itinerant character of the ferromagnetism, pre- 
sumably due to 5/ electrons of uranium. Ulr crystallizes in 
the monoclinic structure (space group P2i) which lacks in- 
version symmetry, and the FM moment is Ising-like, oriented 
along the [101] direction in the (ac)-plane. 

Given the proximity of the SC state in Ulr to the PM tran- 
sition, one may probably consider the magnetization as a 
perturbation on top of the SOC-split bands. Neglecting the 
effect of the former, it is known'^ that even in the case of 
non-centrosymmetric superconductors (and h = 0), the band 
energies still satisfy the relation £^(k) = ep{—'k) due to the 
time reversal symmetry of the single-electron Hamiltonian. 
As a consequence, the SC order parameter on the /3-th sheet 
of the Fermi surface transforms according to one of the irre- 
ducible representations of the normal state point group. In the 
case of Ulr, the point group C2 has two one-dimensional irre- 
ducible representations, denoted A and B. Then the SC order 
parameter is an odd function A(— k) = — A(k) given by2i 
A;^'"(k) (X t(k)(?!>^'"(k), where t(k) is an odd phase fac- 
tor^^ and the basis functions </)^'^ are even in k. Denoting the 
rotation axis of the C2 group as z (this actually corresponds to 
6-axis in case of Ulr), the even functions (j>A and (/j^ can then 
be cast in the following form: 0a (k) = {k1 + C) ui(k), and 
0B(k) = kz{kxU2{k) + kyU^iyCf), where C is some con- 
stant and {Mi(k)} are arbitrary even functions of k^, ky, k^- 
Function (pA generically has no nodes, whereas (j)B has two 
point nodes at the poles (kx=ky=Q) and a line of nodes at the 
equator. The symmetry argument does not allow one to de- 
termine which pairing channel is realized, however the exper- 
imental observation of the strong pair-breaking effect due to 
disorder'^ indicates that the gap must be anisotropic, possibly 
favouring the gap with the nodes such as AbO^)- 

One way of experimentally probing the symmetry of the 
superconducting order parameter in Ulr would be by means 
of transport properties such as Josephson tunneling or point- 
contact spectroscopy. In particular, it has recently been shown 
that the presence of multiple gaps in superconductors with 
broken inversion symmetry should manifest itself through 
clear signatures at bias voltages corresponding to the sum and 
difference of the singlet and triplet components22i2^i2^. We ex- 
pect similar behavior in the present case, at least when the 
ferromagnetism is weak, and point-contact spectroscopy data 
could then be compared with the predictions for Ba in Fig. 
[T^. Alternatively, it should be possible to directly probe the 
spin-texture of the superconducting order parameter by study- 
ing the effect of an externally applied magnetic field when the 
paramagnetic limitation dominates, e.g. in a thin-film struc- 
ture, where the orbital mechanism of destroying superconduc- 
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tivity is suppressed . 

In summary, we have developed a mean-field model for 
a superconductor lacking inversion symmetry and displaying 
itinerant ferromagnetism. Specifically, we have investigated 
the interplay between ferromagnetism and asymmetric spin- 
orbit coupling and how these affect superconducting order, 
which in general is a mixture of a singlet and triplet compo- 
nents. Our main results are the analytical expression Eqs. Q 
and (|4| and the belonging discussion. We find that spin-orbit 
coupling may enhance superconductivity in both the singlet 
and triplet channels in addition to favoring the Stoner crite- 



rion for the ferromagnetic instability. We have applied these 
considerations to the heavy fermion superconductor Ulr, to- 
gether with group-theoretical analysis of the symmetry of the 
SC order parameter. 
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